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optical depth of cloudy media from polarimetric 
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We present the development of a semi-analytical algorithm for optical particle sizing in disperse media. 


The algorithm is applied to the specific case of water clouds. 
modifications to other types of light-scattering medium. 
medium is large and the probability of photon absorption ß is small. 
problem is studied in the regime of highly developed multiple light scattering. 


However, it can be extended with minor 
It is assumed that the optical thickness 7 of the 
Thus the optical particle-sizing 
It was found that the 


degree of polarization in visible and near-infrared channels provides us with information both on the 
effective size of droplets and on the optical thickness t. © 2002 Optical Society of America 
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1. Introduction 


The influence of clouds on radiative fluxes in the 
terrestrial atmosphere has been studied by many au- 
thors both experimentally and theoretically..3= The 
knowledge obtained was applied to the solution of the 
inverse problem, namely, to the determination of the 
thermodynamic phase, microstructure parameters 
(e.g., the effective radius of droplets), single- 
scattering albedo, optical thickness, and liquid water 
path of cloudy media.4—1° 

However, cloudy media do not only absorb and 
scatter radiation. They can also be considered as 
major sources (along with aerosols, gases, and under- 
lying surfaces) of polarized light in the terrestrial 
atmosphere. Clouds introduce a new property of in- 
coming solar light, namely, the preferential direction 
of oscillation of the electric vector in a light beam. 
Generally speaking, both reflected and transmitted 
light becomes elliptically polarized. 

The polarization arises mainly because of the sin- 
gle scattering of a light beam by a water droplet, an 
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aerosol particle, or an ice crystal. Molecular scatter- 
ing is another source of polarized light in the Earth’s 
atmosphere. The phenomenon of multiple scatter- 
ing, which is of particular importance for clouds, 
leads to the randomization of both photon directions 
and polarization states. This causes a decrease of 
light polarization because of an increase in the optical 
thickness ofa cloud. Another effect of multiple scat- 
tering is the introduction of ellipticity in scattered 
light beams. It should be pointed out that the solar 
light scattered by a spherical water droplet becomes 
partially linearly polarized. This is due to the fact 
that the ellipticity of singly-scattered light is equal to 
zero for spherical scatterers if the incident light is 
unpolarized.!! 

Studies of polarization characteristics of solar light 
transmitted and reflected by cloudy media have a 
long history.12-17 However, the real burst of re- 
search in this area was given by a launch of the 
POLDER (polarization and directionality of the 
Earth’s reflectances) instrument on board the Japa- 
nese Advanced Earth Observing Satellite I.18 The 
POLDER was able to transmit to the Earth a huge 
amount of information about polarization character- 
istics of light reflected from cloudy media, aerosols, 
and underlying surfaces at several wavelengths. 
Specifically, the first three components of the Stokes 
vector S,(7, Q, U, V) have been measured for wave- 
lengths à equal to 443,670, and 865 nm. There is no 
doubt that even more advanced polarimeters with 
wider spectral coverage will appear on board differ- 
ent satellites in the future, which makes further the- 


oretical studies of polarization characteristics of 
cloudy media extremely important. This is due to 
potential possibilities for global retrievals of the cloud 
microstructure,’ the shape of particles,!8 and the op- 
tical thickness of clouds?’ on the basis of polarization 
measurements. 

Of course, similar information can be obtained from 
the reflected intensity measurements. However, it 
could well appear that the degree of polarization can be 
used as a source of additional information about the 
cloud particle size distributions close to the top of a 
cloud.2 Indeed, the high proportion of photons scat- 
tered from a thin upper layer of a cloud in the creation 
of light polarization is quite understandable. Multi- 
ply scattered light fluxes from deep layers are hardly 
polarized at all. Radiative characteristics, in con- 
trast, represent the cloud asa whole. Thus the effec- 
tive radius derived from radiative measurements is an 
average of large ensembles, of possibly quite different 
particle size distributions, presented in different parts 
of cloudy media. 

It was emphasized by Hansen!2 30 years ago that 
the fact that the polarization data refer to cloud tops 
may be not a drawback but an advantage. Indeed, it 
is easier to interpret observations that do not refer to 
some average values over the entire cloud depth. 

The polarization characteristics of cloudy media 
can be studied by application of numerical codes2-12 
based on the vector radiative transfer equation solu- 
tion. However, one can also use the fact that cloud 
fields are optically thick in most cases. This allows 
the application of asymptotic analytical relations,!® 
derived for optically thick disperse media with arbi- 
trary phase functions and absorption. These solu- 
tions help us to clear physical mechanisms and main 
features behind the polarization change due to the 
increase of the size of droplets or the thickness of a 
cloud. Analytical solutions also provide an impor- 
tant tool for the simplification of the inverse prob- 
lem.2°. They can be used, e.g., in studies of the 
information content of polarimetric measurements. 

The drawback of the well-known vector asymptotic 
theory for optically thick layers!9 is related to the fact 
that the main equations of this theory depend on 
different auxiliary functions and parameters. Cal- 
culations of these functions and parameters are 
rather complex. This is due to their quite general 
applicability in terms of the single-scattering albedo. 

The special case of weakly absorbing light- 
scattering media was considered by Kokhanovsky?! 
in the framework of the modified asymptotic theory. 
Formulas obtained are rather simple and allow for 
the semi-analytical solution of the inverse problem. 
They can be applied to the case of water clouds. This 
is due to weak absorption of light by water droplets in 
a broad spectral range (V2.2 pm). 

The main task of this paper is to apply equations, 
obtained in the framework of the modified asymptotic 
theory in Ref. 21, to the inverse problem solution, 
namely, to the determination of the droplet size and 
optical thickness of cloudy media. However, results 
obtained after minor changes are made can be ap- 


plied to a much broader range of multiply scattering 
media. 


2. Reflection and Transmission Matrices 


Let us introduce the reflection R(u, wọ, Y) and trans- 
mission Î(u, mo, Y) matrices of cloudy media. Re- 
flection and transmission matrices depend on the 
optical thickness + of a cloud, microstructure of a 
cloud, the solar angle ðo, the viewing angle ð, and the 
relative azimuth | = 6 — po, where $y and ọ are 
azimuths of the Sun and the observer, respectively. 
We also define po = cos 39, p = |cos S|. The impor- 
tance of these matrices is due to the fact that they can 
be used for the calculation of the Stokes vector of 
reflected S,(7, Q, U, V) and transmitted S., Q', U’, 
V') light under arbitrary illumination of a cloudy 
layer. Namely, it follows for plane-parallel cloudy 
layers): 


1 al Qa 
S.(u, b) == awe | do’R(p, p's > 


0 0 


= 6) So(w', $’), (1) 


T pu Qu 
Siu, b) = | dy.’ w’ | dob'T(n, w’, b — ’) 


0 0 
X So(w’, 6’) + Silm, byexo| -*) , (2) 
0 


where So(y’, p’) represents the Stokes vector of the 
incident light beam. The second term in Eq. (2) is 
small for most cloudy media (t >> 1) and can be 
neglected. 

In a good approximation, the Sun is a source of 
light incident from only one direction uo, bo. Then it 
follows: 


1 
Solu, $) = me òlm — Molh — bo) TF, (3) 
0 


where a vector F is normalized so that the first ele- 
ment of nF (namely, mF) is the net flux of solar beam 
per unit area of a cloud layer. Equations (1) and (2) 
simplify if one accounts for Eq. (8): 


Su, o) = Rp, Mo, cu) = bo)F, (4) 
S(p, b) = Tip, Ko, co) — Po)F, (5) 
where only the first element of the column vector 
F: 
_| Pe 
F= F, (6) 
F, 


differs from zero for unpolarized solar radiation. 
Thus only elements R}, T(J = 1, 2, 3, 4) are relevant 
to studies of solar light propagation, scattering, and 
polarization in terrestrial atmosphere. Other ele- 
ments are of interest only in the case of polarized 
(e.g., laser) incident beams. 
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One can see that the knowledge of matrices (u, 
Ho $ — bo) and Êu, mo, $ — oo) is of primary im- 
portance for understanding radiative and polariza- 
tion characteristics of light fluxes that emerge from a 
cloud body. As a matter of fact, they are Mueller 
matrices22.23 for turbid media. Note that Eq. (4) can 
be written in the scalar form 

I= RF, Q= RAF, 


U=R3Fiy, V= Rak, 


(7) 


for the unpolarized incident light. Thus the inten- 
sity of the reflected light is determined by the element 
R,, of the reflection matrix. This is the only element 
that has been a major concern of most cloud optics 


studies up to this date. The degree of polarization is 
given by 

P= JPP +P, (8) 
where 


P,= Ra” T Rg /Ru (9) 


is the degree of linear polarization and 


P, = Ry/Ry (10) 
is the degree of circular polarization. One can 
show?‘ that for particular viewing geometries (e.g., 
uo = 1, p = 1, b — bo = 0, T) the values of Rz; and R41 
are equal to zero and light is partially linearly polar- 
ized. The degree of circular polarization P, is equal 
to zero in this case. It is also common to use a dif- 
ferent definition of the degree of polarization, namely, 


P = -(R/D (11) 


for linearly polarized light beams. It follows from 
Eq. (11) at the nadir-viewing geometry (R31 = R41 = 
0) that 


P= —(Ry/R,,). (12) 


We will use this definition of the degree of polariza- 
tion in this study. The second component of the 
Stokes vector Q is proportional to the difference J, — 
I,, where values of I, and J,, mean the intensities of 
light polarized parallel and perpendicular to a merid- 
ional plane that contains the normal to a layer and 
viewing direction. Thus light polarized perpendicu- 
lar to the meridional plane is characterized by posi- 
tive polarization. Note that the sign of the degree of 
polarization strongly depends on the observation ge- 
ometry. For instance, it is positive in the main rain- 
bow region (see Fig. 1). It could be both positive and 
negative in the glory-scattering region (see Fig. 1). 
Results presented in Fig. 1 were obtained by use of 
the Mie theory for polydispersions of water droplets 
at the wavelength \ equal to 443 nm. The particle 
size distribution was given by the following equation: 


f(a) = Aa’ exp[—(v + 3)a/aeq], (13) 
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Fig. 1. Degree of polarization of light singly scattered by a poly- 
dispersion of spherical water droplets characterized by the gamma 
particle size distribution [Eq. (13)] with the half-width parameter 
v = 6 and the effective radius a.f = 4, 6, 8, 10, and 16 pm as the 
function of the scattering angle at X = 443 nm. (b) Same as in (a) 
but at a.¢ = 6 pm and values of the half-width parameter v = 2, 4, 
6, 8, 10, and 12. 


where a,, is the effective radius of droplets, v is the 
half-width parameter, and 


pr 


A= TG tv)’ 


Note that the mode radius a, is related to the effec- 
tive radius with the following equation: 


3 
av= aul (142), 
Vv 


It should be pointed out that the Stokes vector 
elements J, Q, U, V [see Eq. (7)] can be used for the 
determination of polarization ellipse parameters 
such as the ellipticity and the angle of preferential 


oscillations of the electric vector in a light beam un- 
der study.22 

Thus reflection and transmission matrices allow 
the determination of both the intensity and the po- 
larization states of reflected and transmitted light 
beams under arbitrary illumination of a cloud layer. 


3. Modified Asymptotic Theory 


A. General Equations 


The optical thickness of water clouds is large in most 
cases.! Some mixed and crystalline clouds are also 
extended in the vertical direction for many hundreds 
of meters and could be optically thick.25 Thus the 
asymptotic vector theory!%-26 is quite applicable to 
cloudy media. The reflection matrix R(u, uo, Y, 7) 
and transmission matrix Î(u, uo, T) of a homogeneous 
plane-parallel weakly absorbing turbid layer of a 
large optical thickness T = o,,,L(o.,; is the extinction 
coefficient and L is the geometrical thickness of a 
layer) can be written in the following form?!: 


Rip, Ho, Y, T) = R.(p, Ho, Y) ~~ caine bo, T) 


xX exp(—x — y), (14) 
T(p, Mo, T) = tKo(m)Ko” (uo), (15) 
where 
x=kt, y=4qk, q=[BB0-9]*, 
k =[8(1 — g)(1 - o)l”, 
w=1-B B=, 
O ext 
TS OextL, (16) 


R(u, Wo, Y) is the reflection matrix for a semi-infinite 
medium with the same local optical properties as a 
finite one, L is the geometrical thickness of the me- 
dium, Sext is the extinction coefficient, o,,, is the 
absorption coefficient, g = ; J> p(8)sin 8 cos 6d6 is the 
asymmetry parameter of the phase function p(0), and 
@ is the scattering angle. Note that the vector 
KĮ (uo) in Eq. (15) is transposed to the vector Ko(10). 
The vector K,(j) also occurs in the so-called Milne 
problem.? It describes the angular dependence and 
the polarization characteristics of light emerging 
from semi-infinite turbid layers with sources of radi- 
ation located deep inside the medium. The value of 


TE (17) 
~ sinh(x + ay) 


in Eq. (15) is the global transmittance,?! defined as? 


1 1 
t=4 l udu | PodpoTii(h, Mo, T). (18) 


0 0 


Equation (18) also follows from Eq. (15), taking into 
account the integral relation3 


1 
2 | pdpKo(p) = 1, (19) 


0 


where K,,(w) is the first component of the escape 
vector K,(u). The value of a is approximately equal 
to 1.07.24 

Equations (14) and (15), derived in Ref. 21 for small 
probabilities of photon absorption, are much simpler 
than initial asymptotic formulas, obtained by 
Domke.19 They allow us to calculate the reflection 
and transmission matrices of thick weakly absorbing 
disperse media by a simple means if the solution of 
the problem for a semi-infinite medium with the 
same phase matrix as for a finite layer under consid- 
eration is available.2+?° This reduction of a general 
problem to the case of a semi-infinite medium is of 
general importance for the radiative transfer theory. 


B. Nonabsorbing Media 


It follows from Eqs. (12), (14), and (15) for the degree 
of polarization of the reflected light at the nadir ob- 
servation in the case of nonabsorbing media2!: 


P (ya) 
P = —— 2 
where 
Ko1(1) Koil mo) 
Ulmo) =— a (21) 
ma Ra1(Mo, 1) 
and we have instead of Eq. (17) at B = 0: 
1 
t= i 22 
a + 0.757(1 — g) (aa 
The value of 
R21 (Mo, 1 
P(w) _ 21 (Ko, 1) (23) 


Rai (Mo, 1) 


is the degree of polarization for a semi-infinite non- 
absorbing medium at the nadir observation. Here 
Roy1°(Uo, 1) and R..21°(9, 1) are correspondent ele- 
ments of the reflection matrix [see Eq. (4)] of a non- 
absorbing semi-infinite medium. Note that we 
derived Eq. (20) by taking into account the equality26 
Koo(1) = 0. Owing to this equality, the second com- 
ponent of the Stokes vector for optically thick media 
coincides with the value of Q for semi-infinite layers 
at the nadir observation geometry. This fact can be 
used for the retrieval of the size of particles from 
spectral measurements of Q(\) even if the precise 
information on the optical thickness is not available. 
Thus an accurate estimation of the effective radius of 
droplets can even be obtained from measurements of 
Q at a single wavelength in the visible. 

Both functions Ko; (19) (see Fig. 2) and R,.11°(\1o, 1) 
(see Zege et al.2’, Fig. 3.3b) only slightly depend on 
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Fig. 2. Dependence of the escape function Ko: on the cosine of the 
incidence angle for the Henyey—Greenstein phase function with 
the asymmetry parameter g = 0, 0.5, 0.75, 0.8.28 Results of cal- 
culations with Eqs. (24) are presented by a dashed line. 


the size of scatterers. They can be represented, e.g., 
by the following approximate formulas?’: 


=—(1+ wi (fo, 1) = —-——-. 
Koi(o) 7 (1 + 249), Ra (Ho, 1) 2(1 + po) 
(24) 
Thus it follows from Eqs. (21) and (24): 
54(1 +2 1+ 
U (po) - ( bo) ( Mo) , (25) 


The dependence U(%,) is presented in Fig. 3. One 
can see that U €E [1.08, 1.33]. Thus the function 
U(ðo) does not change drastically with the angle and 
can be substituted by a constant value (e.g., U = 
1.18). Note that the function U(y,) is multiplied by 
the global transmittance t in Eq. (20), the latter being 
a small number for the thick clouds considered here. 
Thus one can see that small errors in the function 


1.35 f i i ' TT] 


10 20 30 40 50 60 70 80 90 
incidence angle, degrees 

Fig. 3. Dependence U(%,) according to Eq. (25). 

represents the average value of U = 1.18. 


The dashed line 
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U(uo) do not much influence the degree of polariza- 
tion calculated with Eq. (20). 

Thus approximate equations (20), (22), (25), and 
(17) allow us to reduce the complexity of calculations 
of the function P(o). Note that exact calculations 
require the solution of the four coupled integro- 
differential equations for the Stokes vector compo- 
nents of the vector radiative transfer theory,® which 
is a far more complex task than application of Eq. 
(20). 

In more general cases [see Eqs. (14) and (15)] the 
functions K,() and R..°(, wo, W) for a nonabsorbing 
semi-infinite light-scattering medium should be 
found numerically. They depend on the phase ma- 
trix of a random medium in question. However, 
they do not depend on the single-scattering albedo wo 
and the optical thickness t by definition. Another 
interesting feature of these functions is their weak 
dependence on the microstructure parameters (size, 
shape, and chemical composition of particles) of a 
medium under study in the broad range of angular 
parameters p, Wo, W. This is due to the randomiza- 
tion of both polarization states and propagation di- 
rections of photons in highly scattering nonabsorbing 
semi-infinite layers irrespective of the local optical 
properties of the medium in which they propagate. 


C. Weakly Absorbing Media 


Equation (20) is easily generalized to account for the 
absorption of light in a cloudy medium. Namely, it 
follows from Eqs. (12) and (14)21: 


P..*(\o) 
1 — U*(po)t exp(—x — y) 


P(o) = > (26) 


where the global transmittance ¢ is given by Eq. (17) 
and 


Ko, (1) K, 
U*(W0) = = (27) 


Values of P„*(uo) and Rs11”(Mo, 1) represent the de- 
gree of polarization and the reflection function of a 
semi-infinite weakly absorbing cloud at the nadir- 
viewing geometry. They can be expressed via the 
correspondent functions for nonabsorbing media. 
Indeed, it follows for the reflection function R, = 
Ren a: 


R.(p, wo, b) = R(p, wo, bexpl—yU(p, po, Y)], 
(28) 


where 


Ulu, mo, $) = Kolu) Kolumo) /R-™ (uo u, $), (29) 


Ky = Ko, and R,’ is the reflection function of the 
semi-infinite nonabsorbing light-scattering medium. 

Our calculations show that the accuracy of Eq. (28) 
can be improved at y > 1 if we substitute the value of 


U by s = (1 — 0.05y)U. Thus we have for the value 
of R.41*(Wo, 1) in Eq. (27): 
Ran*(m, bo, 4%) = tide (ie Lo, Wexp(—sy). (30) 

Substituting Eqs. (27) and (30) into Eq. (26), we 
obtain 

P..* (Wo) 
1 — U(mo)t exp[—x - y(1 - s)]’ 
where we also used Eq. (21). 

Let us find now the relationship between functions 
P.*(\49) and P.,°(9). For this we will use Eq. (12) 
and the fact that the value of R..9;*(uo) for semi- 
infinite weakly absorbing media does not differ much 


from the same value for the nonabsorbing case.24 
Then we have by using Eqs. (12) and (80) 


P.*(o) = P2(u, bo, U)exp(sy). 


4. Inverse Problem 


P(uo) = (31) 


(32) 


A. Algorithm Description 


Equations (20) and (31) can be used for the semi- 
analytical retrieval of droplet sizes from the two- 
wavelength polarimetric inversion algorithm. Let 
us rewrite these equations here, explicitly showing 
the wavelength dependence: 


PLA) 
1 — Ut(d,)’ 


P..°(Az)exp[s(Az)y (Ag) ] 
1 — Ut(Ag)exp{—x (Az) — YADL — s(A2) |} í 
(34) 


where à; is the wavelength in the visible and \, is the 
wavelength in the infrared. We also accounted for 
Eq. (82). A closer look at these equations shows that 
values P(\,) and P(A.) depend only on the solar ele- 
vation, the optical thickness at each wavelength, the 
droplet size distribution, and the refractive index. It 
is known, however, that the influence of the type of 
the particle size distribution on the integral light- 
scattering characteristics, such as x and y in Eqs. (83) 
and (84), is weak.29 So we will neglect this depen- 
dence and characterize the size of droplets by just one 
number—the effective radius of droplets. We also 
ignore the possible variability of the half-width of the 
droplet size distribution and use only the distribution 
[Eq. (13)] at v = 6 in the discussion, which follows. 
The dependence of the degree of polarization of the 
singly scattered light on the value of v is presented in 
Fig. 1(b). We see that this dependence can be ne- 
glected for most scattering angles. It is significant, 
however, at 0 ~ 150—155 degrees. 

The parametrizations for the integral light- 
scattering characteristics [Eqs. (16)] as well as func- 
tions P.P), P(o) at wavelengths \, = 0.65 um 
and A, = 1.55 pm are presented in Appendices A and 
B. We will use these two wavelengths in the devel- 
opment of the inversion algorithm. 


Piha) = (33) 


P(X) = 


degree of polarization, % 


optical thickness 


Fig. 4. Dependence of the degree of polarization on the optical 
thickness at à = 0.65 pm, 1.55 um, So = 37°, Ò = 0°, aef = 6 pm, 
and v = 6 [see Eq. (13)] according to the exact radiative transfer 
calculations (symbols) and Eqs. (33), (34), and (B1). 


Thus we have the system of the two Eqs. (33) and 
(84) to find the optical thickness in the visible and the 
size of droplets. We proceed as follows. First of all, 
we express the value 7(A,) in Eq. (33) via the mea- 
sured degree of polarization P(\,): 


Alt OA) a a] 


= =a ee 
where 
i= 1 Pe'a) (36) 
X U PA) |’? 


and we also used Eq. (22). Second, we use the fact 
that the ratio Y(a.,) = T(\g)/T(A,) depends only on the 
effective radius of droplets (see Appendix A). So we 
can substitute the value of t(Ag) by Yaet) in Eq. 
(34). It gives us a single transcendent equation [Eq. 
(84)] for the effective radius determination. This 
equation can be easily solved with standard tech- 
niques. The value of 7(\,) can be obtained from Eq. 
(85) afterward. This completes the proposed algo- 
rithm description. 


B. Accuracy of the Retrieval Algorithm 


Let us study the errors associated with the retrieval 
algorithm proposed. We will use the fixed viewing 
geometry (ðo = 37°, 3 = 0°), in which the sensitivity 
of the measured degree of polarization in the visible 
to the value of the effective radius of droplets is high 
(see Fig. 1). However, one can also choose another 
geometry. This is mostly due to the fact that the 
information on sizes of droplets is contained in the 
parameter y, which does not depend on the observa- 
tion geometry. This is an important point. 

First, let us study the accuracy of Eqs. (33) and 
(84). The results of comparisons with the exact ra- 
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effective radius oe 
optical thickness _..---~ 


error of retrieval, % 


error of measurements, % 


Fig. 5. Error of the retrieval as the function of the errors of 
measurements. 


diative transfer code are presented in Fig. 4. One 
can see that errors of approximate formulas at optical 
thicknesses larger than 5 are smaller than possible 
measurement and forward-model assumption errors. 
So approximate equations can be applied to the in- 
verse problem solution. 

Let us consider now the error in the retrieval of a.¢ 
and 7+ that is due to possible measurement errors. 
For this, we have calculated the values of P(\,) and 
P(\g) at aef = 6 pm and rt = 10 with the exact radi- 
ative transfer code for wavelengths 0.65 and 1.55 um. 
Then we assumed that the actual values of P(\,) and 
P(\,) differ from the calculated values by the relative 
positive error A, which was changed from 0% to 10%. 
The results are presented in Fig. 5. We see that the 
errors in the retrieved values of a,,and 7 are smaller 
than 12% in this case. This shows that our algo- 
rithm is not particularly sensitive to the measure- 
ment errors at a,, = 6 wm and t = 10. Clearly, the 
errors of retrieval of the optical thickness increase 
with the optical thickness. This is due to the fact 
that the degree of polarization for extremely thick 
clouds is not dependent on the optical thickness at all 
[see Eqs. (383) and (34) at t = 0]. So it cannot be 
retrieved in principle. This is not the case as far as 
the effective radius of droplets is concerned. It can 
be derived for clouds of arbitrarily large thicknesses. 
However, for large clouds the transmission ¢ is ap- 
proximately equal to zero and the transcendent equa- 
tion [Eq. (84)] transforms into more simple equation 
[Eq. (82)]. 

Solving this single transcendent equation numeri- 
cally, one obtains the value of a.ş which subsequently 
can be used to find the optical thickness and liquid 
water path of clouds, if they are not thick. Clouds 
with optical thicknesses around 15 and larger could 
be considered as semi-infinite for wavelengths larger 
than 1.55 pm. Of course, this is not the case in the 
visible. Knowing the values ofa,-and 7, one can find 
it easy to obtain many other parameters of cloudy 
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media, including their transmittance, reflection func- 
tion, and spherical albedo. 

Spectral measurements at several wavelengths 
can improve the accuracy of the algorithm. How- 
ever, remarkably, as we have discussed before, there 
is even a possibility!’ of finding the size of particles 
and the liquid water path of clouds from measure- 
ments of the reflection function and the degree of 
polarization at the single wavelength in the visible. 
For this, one should make measurements at the rain- 
bow geometry. 


5. Conclusions 


We presented here a system of simple equations that 
can be used for parametrizations of radiative and 
polarization characteristics of cloudy media. The al- 
gorithm of the inverse problem solution was pro- 
posed. It is based on spectral measurements of the 
degree of polarization of reflected light in the visible 
and the infrared at the rainbow geometry. Note 
that the method proposed is not significantly influ- 
enced by the choice of the observation geometry. 
This is not the case, e.g., when only measurements in 
the visible are used.9-17 

The reflection of light from the surface and scatter- 
ing by molecules and aerosols was neglected in this 
study. However, they can easily be taken into ac- 
count. The method described can also be used for 
multiple light-scattering disperse media other than 
water clouds. It is also applicable with slight mod- 
ifications for the retrieval of the spectral dependence 
of the imaginary part of the refractive index of drop- 
lets, which can differ from the refractive index of pure 
water because of various sources of contamination. 


Appendix A: Local Optical Characteristics of 
Cloudy Media 


We use the parametrizations of the extinction coeffi- 
cient o,,,, the asymmetry parameter g, and the ab- 
sorption coefficient o,,,,, of cloudy media: 


abs 
3C 1.1 
= 1+ Al 
Oext Qder | ga" , ( ) 
4mxC $ 
Oabs — N > Palta (A2) 
n=0 
4 
=p > uia (A3) 


n=0 


Here C is the volumetric concentration of droplets, xr 
= 2ma,,/\, x is the imaginary part of the refractive 
index of particles. The parameters p, at the 
1.55-u~m wavelength, where there is an appreciable 
amount of light absorption by water droplets, are 
Po = 1.671, pı = 0.0025, pọ = —2.365 x 10+, p3 = 
2.861 x 10 ©, and p4 = —1.05 x 107°. The param- 
eters q„ depend on the refractive index of particles. 
They are given in Table 1 for water droplets at the 
0.65- and 1.55-um wavelengths. The accuracy of 
Eqs. (A1)-(A3) is better than 3% for the effective 


Table 1. Parameters q, for Different Wavelengths à 


A, pm do qı q2 q3 q4 
0.65 0.1121 0.5118 0.8997 0.0 0.0 
1.55 0.0608 2.465 —32.98 248.94 —636 


Table 2. Parameters in Eq. (A4) at the Wavelengths 0.65 and 1.55 pm 


A, pm a S G o 
0.65 9.1 22.0 0.37 0.17 
1.55 5.5 8.4 =1.0 0.29 


degree of polarization, % 


2 4 6 8 10 12 14 16 18 20 
effective radius, um 


(a) 


degree of polarization, % 


4 6 8 10 12 14 16 18 20 22 
effective radius, um 


(b) 
Fig. 6. Degree of polarization, calculated with Eq. (B1) (line) and 
exact radiative transfer code (symbols) at (a) \ = 0.65 um, (b) 1.55 
pm. Other parameters are B = 0, ðo = 37°, 0 = 0°, def = 6 pm, 
and v = 6. 


radius of droplets in the range of 4—16 um, which is 
most common for water clouds. 


Clearly, the ratio of optical thicknesses Y for se- 
lected wavelengths is given by ratios of values {1 + 
[1.1/(«.,)?/7]} for both wavelengths. 


Appendix B: Degree of Polarization at the 
Rainbow Angle 


We present here the results of the parametrization of 
the degree of polarization at the rainbow geometry. 
It is assumed that the incidence angle is equal to 37 
deg and the observation angle is zero, which means 
the nadir observation. Also, we consider the case of 
a semi-infinite nonabsorbing medium. Then the de- 
gree of polarization depends only on the effective ra- 


dius of droplets at this geometry (see Fig. 1). It can 
be presented by the following equation: 
P2=a s (B1) 


1 + expt + oa’ 


where the parameters a, s, č, o are presented in Table 
2 at the 0.65- and 1.55-um wavelengths. The accu- 
racy of these equations is better than 6% at the ef- 
fective radii in the range of 4—16 um (see Fig. 6). 


This research was supported by Volkswagen Fel- 
lowship, Germany. The author is grateful to T. 
Takashima and K. Masuda for providing their vector 
radiative transfer code. 


References 


1. K. Ya. Kondratyev and V. I. Binenko, Impact of Clouds on 
Radiation and Climate (Gidrometeoizdat, Leningrad, 1984). 

2. K. N. Liou, Radiation and Cloud Processes in the Atmosphere 
(Oxford U. Press, Oxford, England, 1992). 

3. I. N. Minin, Radiative Transfer Theory in Planetary Atmo- 
spheres (Nauka, Moscow, 1988). 

4. M. D. King, “A method for determining the single scattering 
albedo of clouds through observation of the internal scattered 
radiation field,” J. Atmos. Sci. 38, 2031-2044 (1981). 

5. M. D. King, “Determination of the scaled optical thickness of 
clouds from reflected soar radiation measurements,” J. Atmos. 
Sci. 44, 1734—1751 (1987). 

6. M. D. King, Y. J. Kaufman, W. P. Menzel, and D. Tanre, 
“Remote sensing of cloud, aerosol, and water vapour properties 
from the moderate resolution imaging spectrometer (MODIS),” 
IEEE Trans. Geosci. Remote Sens. 30, 2—7 (1992). 

7. T. Nakajima and M. D. King, “Determination of the optical 
thickness and effective particle radius of clouds from reflected 
solar radiation measurements. Part 1. Theory,” J. Atmos. 
Sci. 47, 1878-1893 (1990). 

8. Q. Han, W. B. Rossow, and A. Lacis, “Near global survey of 
effective droplet radii in liquid water clouds using ISCCP 
data,” J. Clim. 7, 465-497 (1994). 

9. F.-M. Breon and P. Goloub, “Cloud droplet effective radius 
from spaceborne polarization measurements,” Geophys. Res. 
Lett. 25, 1879-1883 (1998). 

10. A. Arking and J. D. Childs, “Retrieval of clouds cover param- 
eters from multispectral satellite images,” J. Appl. Meteorol. 
24, 323-333 (1985). 

11. H.C. Van de Hulst, Light Scattering by Small Particles (Wiley, 
New York, 1957). 

12. J. E. Hansen, “Multiple light scattering of polarized light in 
planetary atmospheres. Part II. Sunlight reflected by ter- 
restrial water clouds,” J. Atmos. Sci. 28, 1400—1426 (1971a). 

13. J. E. Hansen, “Circular polarization of sunlight reflected by 
clouds,” J. Atmos. Sci. 28, 1515-1516 (1971b). 


20 June 2002 / Vol. 41, No. 18 / APPLIED OPTICS 3657 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


3658 


J. E. Hansen and J. Hovenier, “Interpretation of the polariza- 
tion of Venus,” J. Atmos. Sci. 31, 1137-1160 (1974). 

J. W. Hovenier, “Multiple scattering of polarized light in plan- 
etary atmospheres,” Astron. Astrophys. 13, 7-29 (1971). 

W. J. J. Knibbe, “Analysis of satellite polarization observations 
of cloudy atmospheres,” Ph.D. thesis (Free University of Am- 
sterdam, Amsterdam, Netherlands, 1997). 

A. A. Kokhanovsky, “The determination of the effective radius 
of drops in water clouds from polarization measurements,” 
Phys. Chem. Earth B 25, 471—474 (2000). 

P.-Y. Deschamps, F.-M. Breon, M. Leroy, A. Podaire, A. Bri- 
caud, J.-C. Buriez, and G. Seze, “The POLDER mission: in- 
strument characteristics and scientific objectives,” IEEE 
Trans. Geosci. Remote Sens. 32, 598—614 (1994). 

H. Domke, “Linear Fredholm integral equations for radiative 
transfer problems in finite plane-parallel media: II. Em- 
bedding in a semi-infinite medium,” Astron. Nachr. 299, 95— 
102 (1978). 

A. A. Kokhanovsky, Light Scattering Media Optics: 
and Solutions (Springer-Praxis, Berlin, 2001). 

A. A. Kokhanovsky, “Reflection and transmission of polarized 
light by optically thick weakly absorbing random media,” J. 
Opt. Soc. Am. 18, 883-887 (2001). 


Problems 


APPLIED OPTICS / Vol. 41, No. 18 / 20 June 2002 


22 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


W. A. Shurcliff, Polarized Light (Harvard U. Press, Cam- 
bridge, Mass., 1962). 

H. H. Tynes, G. W. Kattawar, E. P. Zege, I. L. Katsev, A. S. 
Prikhach, and L. I. Chaikovskaya, “Mueller matrix for multi- 
ple light scattering in turbid media: comparison of two inde- 
pendent methods,” Appl. Opt. 40, 400—412 (2001). 

van de Hulst, Multiple Light Scattering (Academic, New York, 
1980). 

A. Uchiyama, S. Asano, M. Shiobara, and M. Fukabori, 
“Ground-based cirrus observations. I. Observation system 
and results of frontal cirrostratus clouds on June 22 and 30, 
1989,” J. Meteorol. Soc. Jpn. 77, 522-551 (1999). 

W. M. F. Wauben, “Multiple scattering of polarized radiation 
in planetary atmospheres,” Ph.D. thesis (Free University of 
Amsterdam, Enschede, Netherlands, 1992). 

E. P. Zege, A. P. Ivanov, and I. L. Katsev, Image Transfer 
through a Scattering Medium (Springer-Verlag, New York, 
1991). 

E. G. Yanovitskij, Light Scattering in Inhomogeneous Atmo- 
spheres (Springer-Verlag, Berlin, 1997). 

J. E. Hansen and L. D. Travis, “Light scattering in planetary 
atmospheres,” Space Sci. Rev. 16, 527—610 (1974). 


